The dynamics of traded value revisited 
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We conclude from an analysis of high resolution NYSE data that the distribution of the traded 
value fi (or volume) has a finite variance (7; for the very large majority of stocks i, and the distribution 
itself is non-universal across stocks. The Hurst exponent of the same time series displays a crossover 
from weakly to strongly correlated behavior around the time scale of 1 day. The persistence in the 
strongly correlated regime increases with the average trading activity (/;) as Hi = Ho + 7 log (fi), 
which is another sign of non-universal behavior. The existence of such liquidity dependent corre- 
lations is consistent with the empirical observation that ai oc {fi)°', where a is a non-trivial, time 
scale dependent exponent. 
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The recent years have seen a number of important con- 
tributions of physics to various areas, among them to fi- 
nance P, [1] . The apphcation of physical concepts often 
seems well suited to the analysis of financial time series, 
however, it is not without caveats. Often, the theoretical 
background of these methods is deeply rooted in physical 
laws that - naturally - do not apply to stock markets. 
In particular, observations regarding power laws 0] , uni- 
versality ^ , and other empirical regularities [s^ are often 
criticized. We carried out a thorough study of the traded 
value per unit time [1, 0, 9 and have arrived at the result 
that some earlier conclusions have to be modified. Here 
we present an analysis of some new data, which supports 
our earlier findings. 

The paper is organized as follows. Section U intro- 
duces notations. Section HIl shows that the distribution of 
traded volume/value is not universal, and it is not in the 
Levy stable regime as sugg estedbyRef. @. Section Im] 
shows, that traded value displays only weak correlations 
for time scales shorter than one day. On longer horizons 
there is stronger persistence whose degree depends loga- 
rithmically on the liquidity of the stock. Finally, Section 
IIVI surveys the concept of fluctuation scaling, shows how 
it complements the observed liquidity dependence of cor- 
relations, and how those two form a consistent scaling 
theory. 



I. NOTATIONS AND DATA 

For a fixed time window size At, let us denote the total 
traded value of the ith stock at time t by 



n,ti{n)e[t,t+At] 



(1) 



where ti(n) is the time when the nth transaction of the 
ith stock takes place. Tick-by-tick data are denoted by 
Vi{n), this is the value traded in transaction n, calculated 
as the product of the price and the traded volume. 

Since price changes very little from trade to trade while 
variations of trading volume are much faster, the fluctua- 
tions of the traded value fi{t) are basically determined by 
those of traded volume. Price merely acts as a weighting 
factor that enables one to compare different stocks, while 
this also automatically corrects the data for stock splits 
and dividends. The correlation properties and the nor- 
malized distribution are nearly indistinguishable between 
traded volume and traded value. 

This study is based on the complete Trades and Quotes 
database of New York Stock Exchange for the period 
1994- 1995. 

Note that throughout the paper we use 10-base loga- 
rithms. 



II. TRADED VALUE DISTRIBUTIONS 
REVISITED 

In this section, we first revisit the analysis done in Ref. 
9]. That work finds that the cumulative distribution 
function of traded volume for time windows of At = 15 
minutes decays as a power-law with a tail exponent A ~ 
1.7 ± 0.1 for a wide range of stocks. This is the so called 
inverse half cube law, and it can be written as 



At 



(2) 
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where Pa* is the probability density function of the same 
quantity. 

The estimation of tail exponents is often difficult due 
to poor statistics of rare events, large stock-to-stock vari- 
ations and the presence of correlations. For the same 
1994 — 1995 period of data and the same 15 minute time 
window certain stocks have A values significantly higher 
than 1.7 [see Fig. [Uleft)]. The tails of these distribu- 
tions can be fitted by a power law over an order of mag- 
nitude, for the top 3 — 10% of the events. The exponent 



2 



A is around 2.8 for these examples. The question arises: 
Which value (if any) is correct? 

In order to address this question we carried out a sys- 
tematic investigation comprising the 1000 stocks with the 
highest total traded value in the TAQ database. We used 
variants of Hill's method [13, EH to estimate the typical 
tail exponent, see Ref. f?! for details. The results of this 
Section are summarized in Table HI Note that in all cases 
the i7-shaped intraday pattern of trading activity was 
removed. 

Most descendants of Hill's method, including the ones 
applied here, contain a free parameter, namely the frac- 
tion p of top events to be considered to belong to the 
tail of the distribution (see Ref. [ill ^''^^l refs. therein). 
According to Fig. [IJleft) this should be set around 
p « 3 - 10%. 

First, let us follow the methodology of Ref. [9] . In that 
paper, the authors first they deduct the mean from the 
time series by taking fi{t) — (fi), where (•) denotes time 
averaging. Then this series is used to estimate value of 
A by applying Hill's method [l3|- The choice p = 0.03 
provides results in line with Ref. for At —15 min 
time windows one finds A = 1.67±0.20. There are several 
issues with this approach: 

1. p is a parameter that can be chosen arbitrarily. 
With the variation of p the same procedure can 
produce estimates ranging from A = 1.1 ± 0.2 
(p = 0.10) to A = 2.15 ± 0.2 {p = 0.005). 

2. The transformation significantly decreases the esti- 
mates of A, down to the range of Levy stable distri- 
butions (A < 2). Estimates for the untransformed 
data are given in Table U for comparison. 

It is simple to show, that the first issue emerges, i.e. the 
estimates systematically depend on p, when one applies 
Hill's method to a finite sample from a distribution of the 
form 

PAt(/) « (/ + /o)-(^+'\ (3) 

where /o is a non-zero constant. The transformation to 
fi{t) — {fi) does not resolve the problem, but biases the 
estimates further. 

Instead, to correct for these biases one can (i) either 
find the proper constant /o, remove it from the data, and 
apply Hill's estimator afterwards (ii) or apply the esti- 
mator of Fraga-Alves [ll[, which is insensitive to such 
shifts. Both of these estimates were found to be signifi- 
cantly higher [2l|: A > 2, see Table HI The methods are 
described in detail in Ref. [3|- 

The two corrected estimators show a strong tendency 
of increasing A with increasing At. Monte Carlo simu- 
lations on surrogate datasets show that this is beyond 
what could be explained by decreasing sample size. For 
distributions with A < 2 increasing window size should 
result in a convergence to the corresponding Levy distri- 
bution, and the measured A's should be independent of 
At. Only when A > 2 can the measured effective value 



of A systematically increase with At. For the 95% of the 
stocks the increasing tendency is observed and for a win- 
dow size At = 1 day the respective A's are greater than 
2. These are strong indications that the distributions 
are not in the Levy stable regime, and thus the second 
moment exists. 

Note that our calculations assume that the variable 
is asymptotically distributed as ^ and do not prove it. 
Still, the existence of the second moment is guaranteed 
by the absence of convergence to a Levy distribution. 
Consequently, it is possible to define the Hurst exponent 
for Mt). 

Regardless of the absence of the convergence to Levy 
stability there are qualitative similarities in the shape 
of the traded value distributions of various stocks [cf. 
Fig. mieft)]. Nevertheless, the existence of a universal 
distribution can be rejected by a simple test [2^ . 

If the form of the normalized distribution was univer- 
sal, then the ratio of the standard deviation and the mean 
would have to obey ct;/ (fi) — h, where h is a constant 
independent of the stock. Equivalently, a relationship 

^^ OC (4) 

would have to hold with an exponent a = 1, at least on 
average. Even though one finds a monotonic dependence 
between the two quantities [as shown in Fig. [Ijright)], 
the exponent is significantly less than 1. This means, 
that the ratio a/ (/) decreases with growing (/), i.e., the 
normalized distribution of / is narrower for larger stocks, 
so their trading exhibits smaller relative fluctuations. We 
will return to this observation in Section [TVl 



III. NON-UNIVERSALITY OF CORRELATIONS 
IN TRADED VALUE TIME SERIES 

One of the classical tools of both financial analysis and 
physics is the measurement of the correlation properties 
of time series [ij, d, [l^. In particular, scaling methods 
[14i] have a long tradition in the study of physical systems, 
where the Hurst exponent Hi is often calculated. For the 
traded value time series /j^*(t) of stock i this is defined 
as 

af(At) = ([.r(t)-(/,^*(t))]^)«At2«', (5) 

Note that it follows from the results of Sect ion [Hit hat the 
variance on the left hand side exists regardless of stock 
and for any window size At. 

The measurements were carried out for all 2474 stocks 
that were continuously available on the market during 
1994 - 1995 nil. Then we sorted the stocks into 6 
groups according to the order of magnitude of their av- 
erage traded value: < (/) < 10^ 10'' < (/) < 10^ . . . , 
10® < (/), all values in USD/min. Finally we averaged 
af{At) within each group. The obtained scaling plots are 
shown in Fig. [21 
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At 


Hill's A (p = 0.06) 


Ref. f9],_p = 0.03 


Shifted Hill's A 


fo/ if) 


Fraga Alves (p = 0.1) 


1 min 


1.43 ±0.09 


1.45 ±0.10 


2.15 ±0.15 


3.0 


1.98 ±0.25 


5 min 


1.56 ±0.13 


1.55 ±0.15 


2.29 ±0.25 


2.8 


2.04 ±0.25 


15 min 


1.71 ±0.20 


1.67 ±0.20 


2.55 ±0.35 


2.8 


2.1 ±0.3 


60 min 


2.06 ±0.30 


1.90 ±0.25 


2.85 ±0.45 


1.8 


2.1 ±0.4 


120 min 


2.3 ±0.4 


2.0 ±0.3 


3.15 ±0.70 


1.6 


2.1 ±0.4 


390 min 


2.7 ±0.6 


2.1 ±0.5 


3.7 ±0.9 


1.2 


no estimate 



Table I: Median of the tail exponents of traded value calculated by four methods for 1994— 1995. The width of the distributions 
is given with the half distance of the 25% and 75% quantiles. 
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Figure 1: (left) The distribution of traded value in At = 15 min time windows, normalized by the mean. The plot displays 
three example stocks for the period 1994 — 1995. The dashed and solid diagonal lines represent power-law distributions with 
tail exponents A — 1.7 and 2.8, respectively, (right) The standard deviation at of trading activity plotted versus the mean {fi) 
of the same quantity. Both for At = 15 min, and At = 390 min « 1 trading day, there is a power law relationship with slopes 
a!(15 min) — 0.74, and a(390 min) = 0.78. A linear proportionality would support the existence of a universal distribution. 
However, the sublinear scaling relationship suggests that relative fluctuations in trading activity are smaller for stocks with 
higher average liquidity. Note: For better visibility, stocks were binned according to (/) and their logo- was averaged. The 
error bars correspond to the characteristic range within the bins. 



All stocks display a crossover around window sizes of 
At = 60 — 390 min, and there are two sets of Hurst ex- 
ponents: valid below, and above the crossover. 
These characterize the strength of intraday and long time 
correlations, respectively. The behavior on these two 
time scales is very different. 

1. For intraday fluctuations, regardless of stock w 
0.51 — 0.52. This means that intraday fluctuations 
of traded value are nearly uncorrelated. 

2. For long time fluctuations the data are correlated, 
but the strength of correlations depends strongly on 
the liquidity of the stock. As one moves to groups of 
larger (/), the strength of correlations (H^) grows, 
up to H+ « 0.8. 

3. If one shuffles the time series, correlations are de- 
stroyed, and i?shuff = 0.5. 



The same phenomenon can be characterized by directly 
plotting the dependence of on (/), as done in Fig. 
[3] Such a dependence is well described by a logarithmic 
law: 



^rf = ffo^+ 7^ log (/,:), (6) 

where 7" = 0.00 ± 0.01, and 7+ = 0.053 ± 0.01. For the 
shuffled time series 7shuff = 0. 

These results indicate, at least in the case of traded 
value, the absence of universal behavior. Liquidity (or, 
analogously, company size) is a relevant quantity, which 
acts as a continuous parameter of empirical observables, 
in particular the strength of correlations and the dis- 
tribution of /. Related results can be found in Refs. 
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log(At) (min) <f> (USD/min) 



Figure 2: The normalized variance ^log(j^(At) — ^ log At 
for the six groups of companies, with average traded values 
(/) G [0, 10*), (/> e [10^ 10^), ...,(/>£ [10«, . . . ) USD/min, 
increasing from bottom to top. A horizontal line would mean 
the absence of autocorrelations in the data. Instead, one ob- 
serves a crossover phenomenon in the regime At = 60 — 390 
mins, indicated by darker background. Below the crossover 
all stocks show very weakly correlated behavior, H~ ~ 0.5. 
Above the crossover, the strength of correlations, and thus the 
slope corresponding to H'^ — i , increases with the liquidity of 
the stock. The asymptotic values of are indicated in the 
plot. 

IV. FLUCTUATION SCALING 

Fluctuation scaling is a general phenomenon, observed 
in a wide range of complex systems [1, H, HI] : A scaling 
law connects the standard deviation ai and the average 
(fi) of the same quantity. In the case of the trading 
activity of stocks we have already presented this result 
in Section |TT] [cf. Eq. ^]. Now we give a more detailed 
discussion. 

Let us start from our observation that 

a.(At) cx , (7) 

where the scaling variable is (fi), or more appropriately 
the stock i, and At is kept constant [see Fig. fright)]. 
Notice that (Ti{At) is the same as in the definition of the 
Hurst exponent in Eq. (O, where i was constant and At 
was varied. 

In Eq. ([7]) the window size At is a free parameter. This 
scaling law persists for any At, but a strongly depends 
on its value, as shown in Fig. [H For small time windows 
(up to 60 min), a(At) « 0.74, then, after a crossover 
regime, when At > 390 min, there is a logarithmic trend. 
This can be summarized as 

a(At) =a± +7±logAt, (8) 

where refers to the regimes At < 60 min and At > 
390 min. The constants are = 0.74, 7" = 0, and 



Figure 3: Value of the Hurst exponents of traded value for the 
time period 1994— 1995. For short time windows (O, At < 60 
min) all signals are nearly uncorrelated, H~ ~ 0.51 — 0.53. 
The fitted slope is 7" = 0.00 ± 0.01. For larger time win- 
dows (■, At > 390 min) the strength of correlations depends 
logarithmically on the mean trading activity of the stock, 
7+ = 0.053 ± 0.01 for 1994 - 1995. Shuffled data (*) dis- 
play no correlations, thus //ghufi = 0.5, which also implies 
7shuff = 0. Note: Groups of stocks were binned, and their 
Hurst exponents were averaged. The error bars correspond to 
the standard deviations in the bins. 

7+ = 0.052 ±0.01. For shuffled time series, a(At) = 0.74 
regardless of At, i.e., 7shuff = 0. 

A visual comparison of Figs. [2] and 2] reveals that the 
crossover in the behavior of a (At) and H falls into the 
same interval. Moreover, when At < 60 min, both Q;(At) 
and H~{{f)) are constant. For At > 390 min, both 
Q!(At) and H~^{{f)) vary logarithmically with their argu- 
ments (see Figs. iniandS]). 

In order to better understand the connection between 
temporal correlations and fiuctuation scaling, let us re- 
peat here Eqs. (O and ([5]): 

a(At) = +7± log At, 
= Ho± + 7* log (/,). 

Beyond the obvious symmetry of these two logarithmic 
laws, notice that the prefactors are equal: in both equa- 
tions 7~ R:! and 7"*" ~ 0.05. 

It is easy to show Q that none of this is a simple 
coincidence. If both fluctuation scaling and long range 
autocorrelations are present in data, there are only two 
possible ways for their coexistence: 

1. Correlations are homogeneous throughout the sys- 
tem. Hi = Hq, 7 = 0, and a is independent of At. 
This is realized for A < 60 min. For shuffled time 
series correlations are absent, thus such data also 
fall into this category. 
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Figure 4: The dependence of the scaling exponent a on the 
window size At for the years 1994 — 1995. The lighter shaded 
intervals have well-defined Hurst exponents and values of 7, 
the crossover is indicated with a darker background. With- 
out shuffling (■) there are two linear regimes: For shorter 
windows a = 0.74 ± 0.02, the slope is 7" = 7(Af < 60 
min) — 0.00 ± 0.01 (solid line), while for longer windows 
a grows logaritmically, with a slope 7^ — 'y{At > 390 
min) = 0.052 ± 0.01 (dashed line). For shuflled data (O) the 
exponent is independent of window size, a{At) — 0.74 ±0.02. 

2. Both the H{{fi)) and a(At) are logarithmic func- 
tions of their arguments with the same coefficient 
7+. This is realized for A > 390 min. 

In other words the coexistence of the two scaling laws is 
so restrictive, that if the strength of correlations depends 
on (/) at all, then the realized logarithmic dependence is 
the only possible scenario. 

V. CONCLUSIONS 

In this paper, we analyzed the empirical properties of 
trading activity on the New York Stock Exchange. We 



showed that, in contrast to earlier findings, the distribu- 
tion of traded value is not in the Levy stable regime, and 
is not universal. Traded value is nearly uncorrelated on 
an intraday time scale, while on daily or longer scales fluc- 
tuations show strong persistence, whose strength grows 
logarithmically with the liquidity of the stock. This ef- 
fect is in harmony with findings on fluctuation scaling, a 
general scaling framework for complex systems. 



All our results imply, that the notion of universality 
must be used with extreme care in the context of fi- 
nancial markets, where the concepts and the theoretical 
background are radically different from those in physics. 
The liquidity of a stock strongly affects the distribution 
and the correlation structure of its trading activity. This 
dependence is continuous, which means the absence of 
universality classes in trading dynamics. The dynamical 
process responsible for such a dependence is yet to be 
identified. 



Finally, we would like to make two remarks. Firstly, 
in Refs. [0] and it is stated that the so called inverse 
half cubic law is observable not only on the 15 minute 
level but also in the tick-by-tick data. Our analysis dealt 
with data aggregated for 1 minute and more, and we 
showed that the assumption of a power law decay is not 
consistent with thejnverse half cubic law for these cases. 
Secondly, in Ref. a footnote mentions the possibility 
of an exponential cutoff in the distribution. This assump- 
tion would influence the estimators strongly and we did 
not consider this case. We thank the anonymous referee 
for calling our attention to these points. 
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